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Abstract. Many observations of pedestrian dynamics, including various self-organization 
phenomena, have been reproduced successfully by different models. But the empirical 
J>-^ databases for quantitative calibration is still insufficient, e.g. the fundamental diagram as one 

i-C ' of the most important relationships displays non-negligible differences among various studies. 

To improve this situation, experiments in straight corridors and T-junction are performed. Four 
different measurement methods are defined to study their effects on the fundamental diagram. 
■ It is shown that they have minor influences for p < 3.5 irT 2 but only the Voronoi method is able 

\^ , to resolve the fine-structure of the fundamental diagram. This enhanced measurement method 

\Q ' permits to observe the occurrence of boundary-induced phase transition. For corridors of 

different widths we found that the specific flow concept works well for p < 3.5 m 2 . Moreover, 
• \ we illustrate the discrepancies between the fundamental diagrams of a T-junction and a straight 

1 corridor. 
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1. Introduction 

During the last few decades, research on pedestrian and traffic flow became popular and 
attracted a lot of attention [1, 2, 3, 4, 5, 6]. The investigation of pedestrian motion plays 
an important role in guaranteeing the safety of pedestrians in complex buildings or at mass 
events. A large number of models have been developed in the past. Most of them are able to 
reproduce phenomena of pedestrian movement qualitatively. Before using a model to predict 
quantitative results like the total evacuation time, it needs to be calibrated thoroughly and 
quantitatively using empirical data. However, this is still difficult due to a lack of reliable 
experimental data. In addition, the small number of available datasets shows surprisingly 
large differences [7, 8]. 

In recent years, several well-controlled pedestrian experiments [9, 10, 11, 12, 13] and 
field studies [14, 15, 16] have been performed. One of the most important characteristics 
of pedestrian dynamics is the fundamental diagram which states the relationship between 
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pedestrian flow and density. Several researchers, in particular Fruin and Pauls [17], 
Predtechenskii and Milinskii [18], Weidmann [19], Helbing et al. [20] have collected 
information about the relation of occupants density and velocity. But there exists considerable 
disagreement among these data. In the comparison performed in [21] the density p , where 
the velocity approaches zero due to overcrowding, ranges from 3.8 m~ 2 to 10 m~ 2 , while the 
density p c where the flow reaches its maximum ranges from 1.75 m 2 to 7 m 2 . Several 
explanations for these discrepancies have been proposed, including cultural factors [22], 
differences between unidirectional and multidirectional flow [23, 24]. 

For single-file movements it was found that even the measurement method has a large 
influence on the fundamental diagram and could be responsible for the observed deviations 
[21]. For identical facilities (e.g., corridors, stairs, doors) with different width it is usually 
assumed that the fundamental diagrams are unified in a single diagram for specific flow J s . 
The study of Hankin et al. [25] in the London subway shows that above a certain minimum 
of about four feet the maximum flow in subway stations is directly proportional to the width 
of the corridor. Besides, it is still not clear whether or not the fundamental diagrams for other 
scenarios like bottlenecks or T-junctions are the same. 

Facing such questions, a series of well-controlled laboratory experiments are carried 
out. The goals of this study are to improve the database related to pedestrian dynamics, to 
determine the influence of measurement methods on the fundamental diagram, and to check 
whether or not the fundamental diagram for different types of facilities can be unified in a 
single diagram. 

In section 2, the experiment setup will be briefly described. The measurement methods 
used in this paper are introduced and defined in section 3. The main results of the work are in 
section 4. Finally, we make some concluding remarks in section 5. 

2. Experiment Setup 

The experiments were performed in hall 2 of the fairground Diisseldorf (Germany) in 
May 2009. They are part of the Hermes project [26] in which the data resulting from the 
experiments will be used to calibrate and test pedestrian movement models. The experiments 
were conducted with up to 350 participants. They were composed mostly of students and each 
of them was paid 50 € per day. The mean age and height of the participants were 25 + 5.7 
years and 1 .76 + 0.09 m, respectively. The free velocity vo = 1.55 ± 0.18 m/s was obtained by 
measuring 42 participants' free movement. 

Figure 1 shows the sketches of the setups and some snapshots during the experiments. 
Two types of geometries, straight corridor (C) and T-junction (T), were used in the 
experiment. 28 runs (see Table 2) were performed in straight corridors with widths of 1.8 m, 
2.4 m and 3.0 m respectively. 7 runs for the T-junction with corridor width of 2.4 m were 
carried out (see Table 3). To regulate the pedestrian density in the corridor, the width of the 
entrance ^entrance and the exit Z? ex it was changed in each run. For details, see Figure 1, Table 2 
and 3. At the beginning, the participants were held within a waiting area. Equal densities 
for different runs were arranged by partitioning the waiting area and counting the number of 
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people in the parts. Standing at the waiting area, they pass through a 4 m passage into the 
corridor. The passage was used as a buffer to minimize the effect of the entrance. In this way, 
the pedestrian flow in the corridor was nearly homogeneous over its entire width. When a 
pedestrian leaves through an exit, he or she returns to the waiting area for the next run. 

The experiments were recorded by two cameras mounted at the rack of the ceiling of the 
hall. To cover the complete region, the left and the right part of the corridor were recorded by 
the two cameras separately. The pedestrian trajectories were automatically extracted from 
video recordings using the software PeTrack [27]. Finally, the trajectory data from the 
two cameras were corrected manually and combined automatically. The frame rate of the 
trajectory data corresponds to 16 fps. Figure 2 shows the trajectories of the head of each 
pedestrian in two runs of the experiment. From these trajectories, pedestrian characteristics 
including flow, density and velocity are determined. 

3. Measurement Methods 

For vehicular traffic it is well known that different measurement methods lead to different 
fundamental diagrams [28, 29]. The results presented in [21] using pedestrian trajectories of 
single file movement, have also shown how large variations induced by different measurement 
methods could be. In previous studies of pedestrian streams, different measurement methods 
were used limiting the comparability of the data. E.G. Helbing et al. proposed a Gaussian, 
distance-dependent weight function [20] to measure the local density and local velocity. 
Predtechenskii and Milinskii [18] used a dimensionless definition to consider different body 
sizes and Fruin introduced the "Pedestrian Area Module" [17]. All of these definitions have 
their advantages and disadvantages. To analyze all of them goes beyond the scope of this 
paper. To enable a detailed analysis, we study the influence of several measurement methods 
on the fundamental diagram and analyze which methods lead to the smallest fluctuations. 

In this study four measurement methods were used to calculate the basic quantities: flow, 
density and velocity. Some terminologies used here are taken from [21] and [30]. 

3.1. Method A 

Method A calculates mean value of flow and density over time. A reference location x in the 
corridor is taken and studied over a fixed period of time At (as shown in Figure 3). We refer to 
this average by () Af . Using this method we can obtain the pedestrian flow J and the velocity 
V; of each pedestrian passing x directly. Thus, the flow over time (J)\ t and the time mean 
velocity (v) Af can be calculated as 



where is the number of persons passing the location x during the time interval At. t NAt 
is the time between the first and the last of the N At pedestrians. Thus t NAt is the actual time 
that the iV Af pedestrians used for passing the location. It can be different from At. The time 
mean velocity (v) Af is defined as the mean value of the instantaneous velocities v,-(/) of the N& t 




and 




(1) 
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(a) Sketch of experiment at straight corridor 



(b) Snapshots from experiment C- 180- 180-070 




(c) Sketch of experiment at T-junction 




(d) Snapshots from experiment T-240- 120-240 



Figure 1. Sketch of the experimental setup ((a) and (c)) and snapshots of the experiment from 
the two cameras respectively ((b) and (d)). The shaded regions in (a) and (c) are the chosen 
measurement areas with 2 m in length. 
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Figure 2. The trajectories for all pedestrians in the experiment. The data are extracted from 
the video recording by the free software PeTrack [27]. 
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Figure 3. Illustration of different measurement methods. Method A is a kind of local 
measurement at cross-section with position x averaged over a time interval At, while Methods 
B — D measure at a certain time and average the results over space Ax. Note that for Method D, 
the Voronoi diagrams are generated according to the spatial distributions of pedestrians frame 
by frame. 
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persons according to equation (2). We calculate v,-(f) by use of the displacement of pedestrian 

i in a small time interval At' around t: 

_ *,•(* + Af72)-s,-('-Af72)) ... 
vKO = ^ (2) 

3.2. Method B 

The second method measures the mean value of velocity and density over space and time. 
The spatial mean velocity and density are calculated by taking a segment Ax in the corridor 
as the measurement area. The velocity (v), of each person is defined as the length Ax of the 
measurement area divided by the time he or she needs to cross the area (see equation (3)), 

Ax 

(v)i = — — (3) 

'out 'in 

where t m and t out are the times a person enters and exits the measurement area, respectively. 
The density p, for each person is calculated with equation (4): 

1 r'° m N'(t) 

{p)i = -t—r • r^k dt (4) 

'out 'in *Jt m u cor 

b cor is the width of the measurement area while N'(t) is the number of person in this area at a 
time t. 

3.3. Method C 

The third measurement method is the classical method. The density <p)a.v is defined as the 
number of pedestrians divided by the area of the measurement section: 

N 

V>** = u T (5) 

Ocor • AX 

The spatial mean velocity is the average of the instantaneous velocities V0) for all pedestrians 
in the measurement area at time t: 

1 N 

<v>a, = -X/<« ( 6 ) 

i=\ 

3.4. Method D 

This method is based on the use of Voronoi diagrams [31] which are a special kind of 
decomposition of a metric space determined by distances to a specified discrete set of objects 
in the space. At any time the positions of the pedestrians can be represented as a set of points, 
from which the Voronoi diagram (see Figure 3) can be generated. The Voronoi cell area, A,-, 
for each person i can be obtained. Then, the density and velocity distribution of the space p xy 
and Vxy can be defined as 

Pxv = l/Ai and v XJ = v,(0 if (x,y) e A t (7) 
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Figure 4. Time series of density and velocity using Method C and D. The two vertical lines 
indicate the beginning and the end of stationary states. 



where v,-(f) is the instantaneous velocity of each person, see equation (2). The Voronoi density 
and velocity for the measurement area is defined as [30] 

<P>, - 4^ (» 
b C0I ■ Ax 

<V>v = J { ~ T (9) 
b C or ■ Ax 



4. Results 



We calculate the fundamental diagram for the straight corridor experiments using the methods 
introduced in the last section. To facilitate a comparison among these four methods, we use 
the hydrodynamic flow equation J = pvb. 



4.1. Influence of the measurement method 

For Method A we choose the time interval At = 10 s, At' = 0.625 s (corresponding to 10 
frames) and the measurement position at x = (see Figure 2). For the other three methods 
a rectangle with length 2 m from x = -2 m to x = and width of the corridor is chosen as 
the measurement area. We calculate the densities and velocities every frame with a frame rate 
of 16 fps. All data below are obtained from the trajectories. To determine the fundamental 
diagram only data at stationary states, which were selected manually by analyzing the time 
series density and velocity (see Figure 4), were used. For Method D we use one frame per 
second to decrease the number of data points and to represent the data more clearly. 

Figure 5 shows the relationship between density and velocity obtained from different 
methods. Using Method A the flow and mean velocity can be obtained directly. To get the 
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Figure 5. The fundamental diagrams, the relationship between density and velocity, measured 
on the same set of trajectories but with different methods. Except for the density in (a), which 
is calculated using p = J/(b ■ Ax), all data are determined directly from the trajectories. The 
legends in (b), (c) and (d) are the same as in (a). 
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relationship between density and flow, the equation p = (J)A t /({v) At ■ b cor ) was adopted to 
calculate the density. For the Method B, C and D the mean density and velocity can be 
obtained directly since they are mean values over space. There exists a similar trend of the 
fundamental diagram obtained using the different methods. However, their influence on the 
scatter of the results is obvious. Table 1 shows the standard deviation of velocities in certain 
density intervals for different measurement methods. Compared to the other approaches, the 
fluctuations in Method B and C are larger. 

Another criterion for the quality of the methods is the resolution in time and space. Even 
if Method A provides a smaller standard deviation than Method D, the low resolution in time 
smears the transition at p = 2m~ 2 clearly visible in Figure 5(d). The density in Method C has 
a strong dependence on the size of the measurement area A m = b C0I ■ Ax. The interval between 
two density values is 1 /A m , which indicates that the measurement area should not be too small 
using this method. But large areas can not provide a detailed spatial resolution. Method D 
can reduce the density and velocity scatter [30]. The reduced fluctuation of Method D is 
combined with a good resolution in time and space, which reveal a phenomenon that is not 
observable with Method A, B and C. In Figure 5(d) it seems that there is a discontinuity of 
the fundamental diagram when the density is about 2 m~ 2 . This will be analyzed in detail in 
section 4.3. 

Figure 6 shows the relationship between the density and flow obtained from different 
methods. The pedestrian flow shows small fluctuations at low densities and high fluctuations 
at high densities. The fluctuations for Method A and Method D are smaller than that for other 
methods. However, there is a major difference between the results. While the fundamental 
diagrams obtained using Method A and Method C are smooth, fundamental diagrams 
obtained with Method B and Method D show a clear discontinuity at a density of about 2 m~ 2 . 
The average over a time interval of Method A and the large scatter of Method C blur this 
discontinuity. In disagreement with the results in [21], no marked differences occur among 
the fundamental diagrams produced by different methods (see Figure 5). In [21], single-file 
movement in a corridor with periodic boundary was studied. In that experiment, distinct stop 
waves occured at high densities and lead to large inhomogeneities in the trajectories. Possibly 
the differences of different methods will be larger in the cases where stop waves occur or when 
the characteristic of the pedestrian flow is not laminar. 



Table 1. Standard deviation of velocities in certain density interval for different methods 



Density interval [m 2 ] 


Method A [m/s] 


Method B [m/s] 


Method C [m/s] 


Method D [m/s] 


pe [0.8,1.2] 


0.119 


0.169 


0.175 


0.120 


pe [1.6,2.0] 


0.086 


0.144 


0.175 


0.111 



4.2. Influence of the corridor width 

From the above analysis it can be concluded that there is no large influence of different 
methods on the results for the density region without the stop wave phenomenon. We have 
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Figure 6. The fundamental diagrams, the relationship between density and flow, measured 
at the same set of trajectories but with different methods. The density in (a) is calculated 
indirectly using p = J/(b ■ Ax), while the flows in (b), (c) and (d) are obtained by adopting the 
equation J = pvb. The legends in (b), (c) and (d) are the same as in (a). 
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shown that the results with the smallest fluctuations are provided by Method A and D. Thus 
we analyze the other unidirectional experiments with corridor width 2.4 m and 3.0 m using 
Method A and Method D. 

Figure 7 shows the relationship between density, velocity and flow using these two 
methods. The fundamental diagram of the same type of corridor but with different widths 
are compared. The fundamental diagrams for these three widths agree well for both methods. 
This result is in conformance with Hankin's findings [25]. He found that above a certain 
minimum of about 4 ft (about 1.22 m) the maximum flow in subways is directly proportional 
to the width of the corridor. Our results agree with the assumption that the specific flow 
J s = J /b is independent of the width of the facility. However, it is possible that for small 
corridors or very high densities J s becomes dependent on b. 

4.3. Discontinuous trend of the fundamental diagram 

In Figure 5(b) and 5(d) a discontinuity occurs at p ~ 2 m~ 2 , separating the function J s (p) 
in a region p < 2 m~ 2 with negative curvature and a region p > 2 mT 2 with positive 
curvature. Moreover, for Method D a gap occurs around v = 0.7 m/s. This transition is 
also found in the experiments with b COI = 2.4 m and b cor = 3.0 m. The fundamental diagram 
changes qualitatively when the width b exit of the exit is modified. The modification of the 
exit width was necessary to achieve high densities. However, it seems that this slight change 
in the experimental setup causes a significant change in the flow-density relation. This point 
becomes obvious in the velocity-density relation, especially in Figure 5(d) which is obtained 
from Method D. Although the measurement area in the corridor is 4 m away from the exit, the 
influence of the change of the exit on the fundamental diagram is sensitive. The decreasing of 
the exit width limits the outflow of pedestrians and leads to a discontinuity in the fundamental 
diagram. This can be interpreted in terms of the well-established theory of boundary induced 
phase transitions, see section 4.4. 

4.4. Interpretation in terms of boundary -induced phase transitions 

Some of the results can be interpreted in terms of the well-established theory of boundary- 
induced phase transitions [6, 32]. In nonequilibrium systems phase transitions (in the bulk) 
can be induced by changing boundary conditions, generically input and output rates in the 
case of transport systems. A mesoscopic theory has been developed which allows to derive 
the phase diagram of an open system (allowing input and output of particles at the boundaries, 
see Figure 8) from the fundamental diagram of the periodic system [33]. This theory even 
makes quantitative predictions on the basis of an extremal principle [34]. 

The phase diagram as function of the boundary rates a and fi has a generic structure. The 
number of phases observed depends only on the number of local maxima in the fundamental 
diagram. For generic traffic systems it has only one maximum and the a-/?-phase diagram 
consists of three phases, the high-density phase (HD), the low-density phase (LD) and the 
maximum current phase (MC), see Figure 9. When the supply rate a of the particles is larger 
than the removal rate /3 and J3 < (3 C , the particle extraction is the limiting process resulting in 
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Figure 7. Comparison of the fundamental diagram in different corridor widths. 



a high density phase where the current is independent of a. When the particles are supplied 
not too fast, a < a c and /? > a, a low density phase is formed which is limited by particle 
supply. Here the current is independent of a. There is a discontinuous phase transition along 
the line a = < a c . When particles are supplied and removed sufficiently rapidly, a > a c 
and (3 > f3 c , a continuous phase transition into a maximum-current phase for which transport 
is bulk dominated [32] occurs where the current is independent of both a and f3. 

The experiments in the corridor geometry can effectively be described by such a scenario. 
The input rate a into the corridor is controlled by the width ^entrance of the entrance whereas 
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Figure 8. Open system with particle input at the left boundary with rate a and particle output 
at the right boundary with rate /3. 




Figure 9. Generic form of the phase diagram for an open system with input rate a and output 
rate /3. In the LD phase the current has the form J = J{a; p) and in the HD phase J = /(/?; p). 
In the MC phase, the current is independent of a and j3 and corresponds to the maximum of 
the fundamental diagram: J = J max (p). 



the output rate j3 is controlled by the width £> ex it of the exit. Therefore the a-jS-phase diagram 
corresponds in our case to a ^entrance -^exit-diagram. The width b cor of the corridor, on the other 
hand, controls the maximal possible bulk flow in the system, given by the maximal flow 7 max 
in the fundamental diagram. 

For the experiments with b ex [ t = b cor the flow through the system is not limited by the 
exit (see Fig. 6), corresponding to the case > fi c . The system is then in the low-density 
phase. Here the flow is controlled by the inflow into the system, i.e. effectively by £ entrance . At 
^entrance ~ 1 .45 m a transition into the maximum current phase can be observed (see especially 
Fig. 6(a)). For ^entrance = ^cor and b ex n < b CO[ the system is in the high-density phase. The 
phase diagram will be investigated in more detail in a separate publication [35]. 
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Figure 10. Comparison of the fundamental diagrams between straight corridor and T-junction 



4.5. Comparison of straight corridor with T-junction 

In this section, we compare the fundamental diagrams for a straight corridor (C) and T- 
junction (T) with channel width b = 2.4 m. 

Figure 10 shows the results obtained from the experiments C and T using Method A. The 
data assigned with T-left' and T-right' are measured in the area before the stream merge (see 
Figure 1(c)). The data assigned with T-front' are measured in the region where the streams 
already merged. The comparison shows that the fundamental diagram of the unidirectional 
flow agrees with the fundamental diagram of the T-junction in the 'front' part (r-front). But 
for the density region p from 0.5 m~ 2 to 2 m~ 2 the velocities at the 'right' and 'left' part of 
the T-junction (r-left and T-right) are significantly lower. For densities higher than 2 m 2 the 
difference is smaller. The differences in the fundamental diagram could be interpreted in the 
following way. The dynamics in the region after the merging of the streams is comparable 
with the unidirectional pedestrian flow in an open corridor. But in front of the merging the 
velocities are significantly lower indicating that the dynamics of the stream changes due to 
the change of the geometry and the merging of the streams. 



5. Summary 

Series of well-controlled laboratory pedestrian experiments were performed in straight 
corridors and T-junction in this study. Up to 350 persons participated in these activities 
and the whole processes of the experiment were recorded using two video cameras. The 
trajectories of each pedestrian are extracted with high accuracy from the video recordings 
automatically using PeTrack. Four measurement methods are adopted in this study and their 
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influences on the fundamental diagram are investigated. It is found that the results obtained 
from different methods agree well and the main differences are the range of the fluctuations 
and the resolution in time of the results. The influence of the corridor width on the results 
is also investigated. It is shown that fundamental diagrams for the same type of facility 
but different widths agree well and can be unified in one diagram for specific flow. From 
the comparison of the fundamental diagrams between straight corridor and T-junction, it is 
indicated that the fundamental diagrams for different facilities are not comparable. The reason 
for this difference may be the equilibration between the inflow and outflow of pedestrians in 
the corridor. When the outflow and the inflow are not equal, a transition between low and 
high densities appears in the pedestrian flow. This transition can also be observed in the 
fundamental diagram. 
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Table 2. Parameters for the straight corridor experiments 



Experiment index 


Name 


^entrance [^hJ 
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